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We have used Monte-Carlo methods and analytical techniques to investigate the influence of the 
characteristic parameters, such as pipe length, diffusion, adsorption, desorption and reaction rate 
constants on the steady-state properties of Single-File Systems with a reaction. We looked at cases 
when all the sites are reactive and when only some of them are reactive. Comparisons between 
Mean-Field predictions and Monte Carlo simulations for the occupancy profiles and reactivity are 
made. Substantial differences between Mean-Field and the simulations are found when rates of 
diffusion are high. Mean-Field results only include Single-File behavior by changing the diffusion 
rate constant, but it effectively allows passing of particles. Reactivity converges to a limit value 
if more reactive sites are added: sites in the middle of the system have little or no effect on the 
kinetics. Occupancy profiles show approximately exponential behavior from the ends to the middle 
of the system. 

PACS numbers: 02.70Uu, 02.60.-x, 05.50.4-q, 07.05.Tp 



I. INTRODUCTION 

Molecular sieves are crystalline materials with open 
framework structures. Of the almost two billion pounds 
of molecular sieves produced in the last decade, 1.4 billion 
pounds were used in detergents, 160 millions pounds as 
catalysts and about 70 millions pounds as adsorbents or 
desiccants. |^ 

Zeolites represent a large fraction of known molecular 
sieves. These are all aluminosilicates with well-defined 
pore structures. In these crystalline materials, the metal 
atoms (classically, silicon or aluminum) are surrounded 
by four oxygen anions to form an approximate tetrahe- 
dron. These tetrahedra then stack in regular arrays such 
that channels and cages are formed. The possible ways 
for the stacking to occur is virtually unlimited, and hun- 
dreds of unique structures are known. |^ 

The channels (or pores) of zeolites generally have cross 
section somewhat larger than a benzene molecule. Some 
zeolites have one-dimensional channels parallel to one an- 
other and no connecting cages large enough for guest 
molecules to cross from one channel to the next. The 
one-dimensional nature leads to extraordinary effects on 
the kinetic properties of these materials. Molecules move 
in a concerted fashion, as they are unable to pass each 
other in the channels. These structures are modeled 
by one-dimensional systems called Single-File Systems 
where particles are not able to pass each other. A par- 
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tide can only move to an adjacent site if that site is not 
occupied. 

This process of Single-File diffusion has different char- 
acteristics than ordinary diffusion which affects the na- 
ture of both transport and conversion by chemical re- 
actions. For Single-File diffusion, the mean-square dis- 
placement of a particular particle is proportional to the 
square-root of time 

(r^) = 2Ft^ 

where F is the Single-File mobility. |^ This is in contrast 
to normal diffusion, where mean-square displacement is 
directly proportional to time. A variety of approaches 
have been used to describe the movement of the particles 
in Single-File Systems, most of them concentrated on the 
role of the Single-File diffusion process. 

Molecular Dynamic(MD) studies of diffusion in zeo- 
lites have become increasingly popular with the advent of 
powerful computers and improved algorithms. In a MD 
simulation the movement is calculated by computing all 
forces exerted upon the individual particles. MD results 
have been found to match experimental observations of 
Single-File diffusion for systems with one type of molecule 
without conversion and with very short pores. ^ |6[ |^ 
Because a molecule can move to the right or to the left 
neighboring site only if this site is free, MD simulations 
under heavy load circumstances require a high computa- 
tional effort for particles that hardly move. However, the 
level of detail provided by MD simulations is not always 
necessary. 

Thus, deterministic models are used also but they 
are mainly focused on dynamic and steady-state infor- 
mation of short pore systems. |^, ^ Several re- 
searchers ^ ^ used a stochastic approach, i.e.. 
Dynamic Monte Carlo(DMC), to determine the proper- 
ties of Single-File Systems. In DMC reactions can be 
included. The rates of the reactions determine the prob- 
ability with which different configurations are generated 
and how fast (at what moment in time) new configura- 
tions are generated. The most severe limitation of the 



DMC method arises when the reaction types in a model 
can be partitioned into 2 classes with vastly different re- 
action rates. In this case, extremly large amounts of com- 
puter time are required to simulate a reasonable number 
of chemical reactions. However, in general the system 
can be simulated for much longer times than with MD. 

All the previous references put the emphasis on the 
transport properties of adsorbed molecules as the impor- 
tant factor in separation and reaction processes that take 
place within zeolites and other shape-selective microp- 
orous catalysts. Rodenbeck and Karger solved numer- 
ically the principal dependence of steady-state properties 
such as concentration profiles and the residence time dis- 
tribution of the particles, on the system parameters for 
sufficiently short pores. In multiple papers, Auerbach et 
al. ^ used Dynamic Monte Carlo to show different 
predictions about Single-File transport and direct mea- 
surements of intercage hopping ion strongly adsorbing 
quest-zeolite systems. Saravanan and Auerbach |l^ 
studied a lattice model of self-diffusion in nanopores, to 
explore the inflvience of loading, temperature and adsor- 
bate coupling on benzene self-diffusion in Na-X and Na-Y 
zeolites. They applied Mean-Field(MF) approximation 
for a wide set of parameters, and derived an analyti- 
cal diffusion theory to calculate diffusion coefficients for 
various loadings at fixed temperature, denoted as "diffu- 
sion isotherms" . They found that diffusion isotherms can 
be segregated into subcritical and supercritical regimes, 
depending upon the system temperature relative to the 
critical temperature of the confined fluid. Supercriti- 
cal systems exhibit three characteristic loading depen- 
dencies of diffusion depending on the degree of degen- 
eracy of the lattice while the subcritical diffusion sys- 
tems are dominated by cluster formation. Coppens and 
Bell ^ studied the influence of occupancy and 

pore network topology on tracer and transport diffusion 
in zeolites. They found that diffusion in zeolites strongly 
depends on the pore network topology and on the types 
and fractions of the different adsorption sites. MF cal- 
culations can quickly estimate the diffusivity, although 
large deviations from the DMC values occur when long- 
time correlations are present at higher occupancies, when 
the site distribution is strongly heterogeneous and the 
connectivity of the network low. 

Few researchers included also reactivity in Single-File 
Systems. Tsikoyannis and Wei |^ considered a reac- 
tive one-dimensional system with all the sites reactive 
in order to get more information about the reactivity 
and selectivity in one-dimensional systems. They used a 
Markov pure jump processes approach to model zeolitic 
diffusion and reaction as a sequence of elementary jump 
events taking place in a finite periodic lattice. Monte 
Carlo and approximate analytical solutions to the de- 
rived Master Equation were developed to examine the 
effect of intracrystalline occupancy on the macroscopic 
diffusional behavior of the system. One conclusion was 
that better results using analytical approach can be ob- 
tained compared to DMC simulation results by including 



more correlations between neighboring sites in regions of 
the systems with high occupancy gradients and less cor- 
relations in regions with low and no occupancy gradients. 
Starting from Wei results about correlations in Single- 
File Systems, Okino and Snurr ]lO| used a deterministic 
model where each site was assumed to have equal activ- 
ity towards reaction. Doublet approximation was found 
to overpredict the occupancy of the sites and the increas- 
ing mobility raised the concentration of reactants in the 
pore. 

Using DMC simulations we have observed that even for 
infinitely fast diffusion, we still have Single-File effects in 
the system. Instead of focusing on diffusion at different 
occupancies of the system, we therefore concentrate in 
this paper on the reactivity of the system, studying the 
reactivity of the system for different sets of kinetic pa- 
rameters, the length of the pipe and the distribution of 
the reactive sites. We analyse the situations when MF 
gives good results and when MF results deviate strongly 
from the DMC simulations. We investigate the effect of 
the various model assumptions made about diffusion, ad- 
sorption/desorption, and reaction on the overall behav- 
ior of the system. We look at the total loading, loading 
with different components, generation of reaction prod- 
ucts and occupancies of individual sites as a function of 
the various parameters of a Single-File System. 

In section II we specify our mathematical model for 
diffusion and reaction in zeolites together with the theo- 
retical background for the analytical and simulation re- 
sults. In section [II A we present the various results for 



the simplified model without conversion. In section III B 
we use MF theory to solve the Master Equation govern- 
ing the system behavior for the case when all the sites 
have the same activity towards conversion. Similarly the 
results obtained using DMC simulations are presented in 
section III B 2 and are compared with MF results. We 
pay special attention to the infinitely fast diffusion case 
and to the influence of the length of the pipe on the over- 
all behavior of the system. In section p^II C we analyze 
again the MF and simulation results but for the case 
when only some of the sites are reactive. The influence 
of the position and number of reactive sites on the reac- 
tivity and site occupancy of the system is outlined. The 
last section summarizes our main conclusions. 



II. 



THEORY 



In this section we will give the theoretical background 
for our analytical and simulation results. First we will 
specify our model and then we will show that the defined 
system obeys a Master Equation. We will simulate 
the system governed by this Master Equation using DMC 
simulations. The rate equations used for the derivation 
of the analytical results are outlined. 



A. The Model 



B. Master Equation 



Because we are interested in reaction of molecules in 
Single-File Systems, we call the system we are modelling, 
Single-File System with conversion. We model a Single- 
File System by a one-dimensional array of sites, each 
possibly occupied by a single adsorbate. The sites are 
numbered 1, 2, . . . , S. An adsorbate can only move if 
an adjacent site is unoccupied. The sites could be reac- 
tive or unreactive and we note with Nreac the number of 
reactive sites. A reactive site is the only place where a 
reaction may take place. 

We consider two types of adsorbates, A and B, in 
our model and we denote with X the site occupation 
of a site, X=(*, A, B), which stands for an empty 
site, a site occupied by A, or a site occupied by a B, 
respectively. The sites at the ends of the system are 
labeled with m, and the reactive sites are labeled with 
r (see figure |l|). We restrict ourselves to the following 
mono and bi-molecular transitions. 

a) Adsorption and desorption 

Adsorption and desorption take place only at the 
two marginal sites i.e., the left and rightmost sites at 
the ends of the system. 



A (gas) -I- *„i — > A„j 
A,„ — > A(gas) + *m 
B„ — > B(gas) + 

where m denotes a marginal site. Note that there is no 
B adsorption. B's are formed only by a reaction. 

b) Diffusion 

In the pipe, particles are allowed to diffuse via hopping 
to vacant nearest neighbor sites. 



Reaction kinetics is described by a stochastic process. 
Every reaction has a microscopic rate constant associated 
with it that is the probability per unit time that the 
reaction occurs. Stochastic models of physical systems 
can be modelled by a Master Equation. j2l] 

By a, /3, we will indicate a particular configuration of 
the system i.e., a particular way to distribute adsorbates 
over all the sites. Pa{t) will indicate the probability of 
finding the system in configuration a at time t and Wap is 
the rate constant of the reaction changing configuration 
[3 to configuration a. 

The probability of the system being in configuration 
a at time t + dt can be expressed as the sum of two 
terms. The first term is the probability to find the sys- 
tem already in configuration a at time t multiplied by 
the probability to stay in this configuration during dt. 
The second term is the probability to find the system in 
some other configuration (3 at time t multiplied by the 
probability to go from P to a during dt. 

P^{t + dt) = (l-dt^W^„)P„(i)+dt^W^„0P^(t) 

P P 

(1) 

By taking the limit dt ^ this equation reduces to a 
Master Equation: 



dPjt) 
dt 



J2[Wa.(3Ppit)-Wf,Mt)]. (2) 

P 



Analytical results can be derived as follow. The value 
of a property X is a weighted average over the values Xa 
which is the value of X in configuration a: 



{X) =J2Pc.Xo.. 



(3) 



A„ + *n+i 

B„ + *„+j 



*n + B„+j , 



where the subscripts are site indices: n—1, 2, ... , S-1. 
c) Reaction 

An A can transform into a B at a reactive site. 



From this follows the rate equation 

d{X) _^dP^^ 
dt ^ ^ dt " 

a 

= Y,[w^pPp -WpMx^ 

= ^VF„;3P^(X,-X^). 

a/3 



(4) 



Br 



The initial state of the system is all that all sites are 
empty (no particles in the pipe). In this paper we will 
only look at steady-state properties and not to the time 
dependence of the system properties starting with no par- 
ticles. 



C. Dynamic Monte Carlo 

Because it might be not always possible to solve the 
Master Equation analytically, DMC methods allow us to 
simulate the system governed by the Master Equation 




FIG. 1: Picture of a Single-File System with two types of adsorbates, A(lighter colored) and B(darker colored). The marginal 
sites are labeled with m, and the reactive sites(lighter colored) with r. Adsoption of A and desorption of A and B can take 
place only at the two marginal sites. An A can transform into a B only on r labeled sites. 



over time. We simplify the notation of the Master Equa- 
tion by defining a matrix W containing the rate constants 
Wapi and a diagonal matrix R by R^js = '^-yW^p, if 
a = [3, and otherwise. If we put the probabilities of 
the configurations Pa in a vector P, we can write the 
Master Equation as 



(R W)P 



(5) 



where R and W are assumed to be time independent. 
We also introduce a new matrix Q, Q(t) = exp[— Rt]. 

This matrix is time dependent by definition and we 
can rewrite the Master Equation in the integral form 



DMC method is not to compute probabilities Pa{t) ex- 
plicitly, but to start with some particular configuration, 
representative for the initial state of the experiment one 
wants to simulate, and then generate a sequence of other 
configurations with the correct probability. The method 
generates a time t' when the first reaction occurs accord- 
ing to the probability distribution 1 — exp[— i?Q,Qt]. At 
time t' a reaction takes place such that a new configura- 
tion a' is generated by picking it out of all possible new 
configurations (3 with a probability proportional to Wa'a- 
At this point we can proceed by repeating the previous 
steps, drawing again a time for a new reaction and a new 
configuration. 



V(t) = Q(i)P(0) -f / dt'Qit - t')WP{t'). (6) 
Jo 

By substitution we get of the right-hand-side for P{t') 

p{t) = m) 

+ [ dt'Q{t-t')WQ{t') 

t t' (7) 
+ f dt' f dt"Cl{t - t')WQ{t' - t")WQ{t") 
Jo Jo 

+ ...]P(0). 

Suppose aX t = the system is in configuration a 
with probability Pa{0). The probability that, at time 
t, the system is still in configuration a is given by 
Qaait)Pa{0) = exp{~ Raat)Pa{0) ■ This shows that the 
first term represents the contribution to the probabilities 
when no reaction takes place up to time t. The ma- 
trix W determines how the probabilities change when 
a reaction takes place. The second term represents the 
contribution to the probabilities when no reaction takes 
place between times and t' , some reaction takes place at 
time t' , and then no reaction takes place between t' and t. 
The subsequent terms represent contributions when two, 
three, four, etc. reactions take place. The idea of the 



III. RESULTS AND DISCUSSION 

A. No conversion 

We mention in this section various results for the sys- 
tem without conversion. These results can be derived 
analytically. The derivations are not difficult, so for com- 
pleteness we give them in the appendix. We will use the 
results when we deal with the system with conversion. 

In a Single-File System without conversion, the rele- 
vant processes to describe are adsorption, desorption and 
diffusion. So, Wap is given by 



where A^'^ equals 1 if a reaction of type "rx" can trans- 
form the system from f3 to a, and equals otherwise. 
Wads, Wdcs, Wdis are the rate constants of adsorption, 
desorption and diffusion respectively. 

If we substitute expression (8) into the Master Equa- 



tion (2), we get 



and 



dt 



(ads) 



+ 



a0 ^'0 



A (ads) p 



Pb-A 



(des) ■ 



73a 

. (diff) 
^/3a 



(9) 



P. 



Using this expression we can show that when the sys- 
tem is in steady state then the prohahiUty of finding the 
system in a certain configuration depends only on the 
number of particles in the system. 



Pa = q{Na) 



(10) 



where Na is the number of particles in configuration a. 
The expression for q{N) is: 



q{N) = 



Wdes 


5 


'Wads' 


N 


Wdos + Wads. 




[WdesJ 





(11) 



Note that diffusion has here no effect on steady-state 
properties. 

The loading of the pipe, defined as the average number 
of particles per site, is then 



1 ^ 



W 



ads 



JV=0 



Wads + Wde 



(12) 



where p{N) is the probability that there are N particles 
in the system. Note again that diffusion doesn't influence 
the steady-state loading. 

The standard deviation, i.e., the fluctuation in the 
number of particles is then: 



(13) 




(Wdes + Wads)2 



To determine how the parameters of the system influ- 
ence the kinetics of the system, we are interested in the 
correlation in the occupancy between neighboring sites. 
We look at one site occupancy and at two sites occupan- 
cies. We denote by (A„) the probability that an A is at 
site n and with (A„A„+i) the probability to have an A 
at site n and one at site n + 1. 

One and two-site probabilities can be derived from the 
fact that all configurations with the same number of par- 
ticles have equal probability and the expressions for q{N). 
We find 



(A„) 



ads 



Wads + Wde 



(14) 



(A„A„+i) = 



Wads 



w. 



ads 



Wdos 



(15) 



Note that this probability does not depend on the site, 
all sites have equal probability to be occupied and that 
there is no correlation between the occupation of neigh- 
boring sites. Again diffusion doesn't influence these prop- 
erties. Note also that these expressions are the same as 
for a model in which particles are allowed to pass each 
other. 



B. All sites reactive 

We look first at the situation with all sites reac- 
tive: i.e., conversion of an A into a B particle can 

take place at any site including the marginal sites. For 
simplicity we consider WdesA=WdosB=Wdos7 and also 
WdiffA=WdiflFB=Wdiff. We will be looking at the total 
loading (Q), the total loading of A's (Qa), the total load- 
ing of B's (Qb), the mrniber of B's produced per unit time 
{Bprod}, and how the distribution of A's and B's varies 
from site to site ((A„) and (B„)). 

Note that the total loading of the pipe for the model 
with conversion is the same as for the model without 
conversion 

Wads 



Q = 



Wads + Wde 



(16) 



The loadings and the production of B's can easily be 
derived from the probabilities (A„) and (B„) so we first 
focus on them. For a non-marginal site we can write 

^=i?„(A,diff)^^(.x)^ (17) 

where i?„'-^''^'*-' is the rate of diffusion of A from and 
to site n, and Rn^'^^^ is the rate of conversion of A to 
B on site n. The conversion takes place at one site and 
is therefore easier to handle than the diffusion. Using 
equation (4) we have 

i?„(-) = W,^ Yl Aa/™^P^(A„c« - A„^), (18) 

where Ana = 1 if site n is occupied by an A in configu- 
ration a and A„q, = if not. We have An,,. — An3 if 
there is an A at site n in configuration /? (A„^=l) that 
has reacted to a B leading to configuration a (A„c = 0). 
This gives us 



r> (rx) 



-Wrx^ Pp = -W,x(A„), (19) 





where the prime restricts the summation to those /3's 



with A 



ri/3 



1. For the diffusion we similarly get 



i?„(A.diff) = w^.^ A„/A''^'ff)P;3(A„„ - An/s). (20) 




There are four ways in which A„ct — A„^ ^ and 
Aq,/?'-^''*'*^'' ^ in /3: there is an A at site n that can move 
to site (n— 1), there is an A at site n that can move to 
(n + 1), there is an A at site in—X) that can move to site 
n and there is an A at site (n+ 1) that can move to site n. 
In all cases we have iS.ap^^'^^'^^ = 1. In the first two cases 
we have A„q, — A„^ = — 1 and in the last two we have 
Ana — A„^ = 1. The summation over /3 in the first case is 
restricted to configurations with an A at site n and a va- 
cant site (n — 1). This gives a term — Wdiff(*n_iA„). The 
other cases give terms -Wdifi(A„*„+i), Wdiff(A„-.i*„), 
and Wdiff (*nA„-|-i). The rate equations then becomes 



dt 



Wdiff'[-(A„*„_i) - (*„_iA„> + (A„_i*„) 



+ (*„A„+i)] - Wrx(A„). 



(21) 



For (B„) we get similarly 
d(B„ 



dt 



= Wdiff'[-(B„*„_i) - (*„_iB^ 
+ (*„B„+i)] + W^rx(A„). 



(B. 



(22) 



The marginal sites have also adsorption and dcsorption. 
They can be dealt with as the conversion. The rate equa- 
tions for A are 



rf(Ai 
dt 



d{As) 
dt 



^ Wdiff [-(Ai*2) + (*lA2)] + Wads(*l) 

- W^diffh(As*S-l) + (*sAs-i)] + VKads(*s) 



M^dcs(As) - W„(As), 



(23) 



and the rate equations for B 
d{B, 



dt 
d{Bs) 
dt 



- I^diff [-(Bl*2> + (*lB2)] - W^dcs(Bl) + M^rx(Ai), 

- W^diff[-(As*s-i) + (*sAs-i)] - I^dcs(Bi) 



(24) 



Note that these coupled sets of differential equations arc 
exact. 



1. Mean Field results 

We will now look at the loadings Qa and Qb and the 
site occupation probabilities (A„) and (B„). We will 
first determine steady-state properties using the (MF) 
approximation: i.e, we put (A„*„_|_i)=:(A„) etc. 
in the rate equations. This gives us 



KB„+i) + (B„_i) - 2(B„)] + W'„(A„), 

[(A2) - (Ai)] - lyrx(Ai) - Wdes(Ai) 



[(A5-1) - (As)] - W„(As) - Wdcs(A5) 



WdisWdcs 
Wads + Wdc 

Wads + Wdc 
_^ WidsWdcs 
Wads + Wdes ' 

„ _ WdiffWdcs 
~ Wads + Wdc 
WadsWdcs 
Wads + Wdcs' 

l^ads + vVdcs 

(25) 

We have used here the probability for a site to be vacant 
that we have determined for the case without conversion. 

We note that these equations are identical to the MF 
equations of a system in which the particles can move 
independently with a rate constant for diffusion equal 
to WdiffWdcs /(Wdes + Wads). This mcaus that the MF 
does not really model the non-passing that characterizes 
a Single-File System. 

The continuum limit of the MF equation is 



(da/dt\ _^fl~b a \ fd^a/dx^ , ^ ... 
[db/dt [ b 1-a) [dH/dx^ ' + 



(26) 



where a=a{x, t) is the probability distribution of A's (a 
similar definition holds for b), and D=Wdisd^ , with d the 
distance between neighboring sites (see appendix). These 
are the equations that are normally used to describe dif- 
fusion in Single-File Systems, [^g || ||, || As this 
equation is derived from the MF equations, it has the 
same drawback; i.e., the Single-File behavior is only in- 
corporated by the reduction of the diffusion, but it does 
effectively allow for passing of particles. This shows up 
as so-called counter diffusion of A's and B's. 



|22, 23 



We see that equations (25) are linear and we can solve 
them at least numerically. We think however that it is 
worthwhile to use an analytical approach. We consider 
the ansatz 

(A„) cx (27) 
in the steady-state equations (25) for (A„). This leads to 
x^ -2{l + a)x + 1^0, (28) 

with 

Wrx Wdes + Wads 



2Wdiff Wdc 



(29) 




FIG. 2: The characteristic length A as a function of xi. 



FIG. 3: Bprod per unit time at one marginal site as a function 
of W^d. for S=Nreac=30, Wdcs=0.8, Wms = 2 and Wrx = 0.4. 



The quadratic equation yields two solutions xi and X2 
with X2 = Xi^ . We have xi — X2 — 1 only when a = 0, 
i.e. when Wrx = 0. We will therefore assume a > and 
xi < 1. Then 



xi = (1 + a) - vMa + 2). 



(30) 



We can write then the solution (A„) = ai(a;i)" + 
a2{xi)^~^^~"' . The symmetry in the occupancy of the 
pipe (A„) — (As+i_,i) yields oi = a2 = a. So, the gen- 
eral solution for the steady state has the form: 



(31) 



The coefficient a is to be determined from the equa- 
tions for the marginal sites in the set of equations(25). 
In the left side of the system n is small and {S +l — n) is 
large. Because xi < 1 we can neglect the second term in 
equation (21) and (A„) cx xi". This means that the prob- 
ability of finding an A at site in the left-hand-side of the 
system is an exponentially decreasing function of the site 
index. If we write (A„) oc e~'s , we find A = —l/ln{xi) 
for the characteristic length of the decrease. The loga- 
rithm makes this length only a slowly varying function 
of the rate constants(see figure ||). When Wdiff becomes 
larger, a approaches 0, xi approaches 1 and A diverges. 
Note that this is a MF result. We will see that in the 
simulations A remains finite. Also when the conversion 
is slow more A's are found away from the marginal sites. 
The second factor in the expression for a equals the re- 
ciprocal of a site being vacant. Low loading leads to a 
smaller a than high loading. Because of the vacancies 
the A's can penetrate farther into the system before be- 
ing converted. For slow conversion or fast diffusion a is 
small and A can be approximated by 



A = 



Wdiff Wde 



The total loading with A's, Qa, is 



(32) 



(33) 



set 


Wads 


Wdcs 


Wdiff 


W.X 


a) 


0.2 


0.8 


0.05 


0.01 


b) 


0.2 


0.8 


0.05 


0.1 


c) 


0.2 


0.8 


2 


0.1 


d) 


0.2 


0.8 


1 


2 


e) 


0.2 


0.8 


10 


2 


f) 


0.8 


0.2 


0.05 


0.01 


g) 


0.8 


0.2 


0.05 


0.1 


h) 


0.8 


0.2 


2 


0.1 


i) 


0.8 


0.2 


1 


2 


j) 


0.8 


0.2 


10 


2 



TABLE I: The sets of parameters used for the simulations 



so, the expression for Qa is 



^ n=l 



2a xi(l - xi"^) 
'J 1-Xi 



Qb = Q - Qa- 
The total production of B's is 

Bprod = WrxQA-S" = 2aWr 



Xi(l - Xi'^) 
1 — Xi 



, (34) 



(35) 



(36) 



2. Simulation results 

We present now the results for different sets of param- 
eters and we compare them with MF results. Because 
we can see from equation (36) that larger pipes don't in- 
crease the productivity of the system, we consider for the 
comparisons of the results a system size S = 30. We have 
considered separately the sets of parameters in Table |. 

The sets of parameters from a) to e) are for the cases 
of low loading and from f) to j) for the high loading. The 
parameters in the table describe the following situations: 
a) and f) for very slow reaction and slow diffusion; b) and 
g) for slow reaction and slow diffusion; c) and h) for slow 
reaction and fast diffusion; d) and i) for fast reaction and 







Qa 




^prod 


Q 


set 


MF 


Sim 


MF 


Sim 


Sim 


a) 


0.0330 


0.0318 


0.0099 


0.0100 


0.209 


b) 


0.0149 


0.0148 


0.0491 


0.0472 


0.198 


c) 


0.0385 


0.0342 


0.1156 


0.1024 


0.204 


d) 


0.0040 


0.0041 


0.2449 


0.2463 


0.200 


e) 


0.0046 


0.0044 


0.2767 


0.2729 


0.201 


f) 


0.0798 


0.0748 


0.0239 


0.0235 


0.795 


g) 


0.0376 


0.0373 


0.1129 


0.1157 


0.804 


h) 


0.0598 


0.0486 


0.1796 


0.1406 


0.802 


i) 


0.0048 


0.0049 


0.2931 


0.2943 


0.797 


j) 


0.0050 


0.0049 


0.3013 


0.2957 


0.801 



For the case W^ds 



TABLE II: Simulation and MF results for Qa and Bp^od for 
all the sets of parameters 



slow diffusion; e) and j) for fast reaction and fast diffu- 
sion. 

We can see from Table || that the simulation and MF 
results match for all the cases except the cases when we 
have low reaction rates and fast diffusion for both low 
and high loading. In these cases MF overestimates the 
amount of A's in the pipe, and consequently overesti- 
mates the B production. In figure || we have the site 
occupancy with A and B both from the simulations and 
MF. We again see that the MF and the simulation re- 
sults agree reasonably well, except for low reaction rates 
and fast diffusion. MF overestimates the characteristic 
length A and allows A's to penetrate farther into the 
pipe than in the simulations. The reason for this is that 
MF describes the fact that the particles cannot pass each 
other by reducing the diffusion, but this effectively does 
allow for passing. The larger A in MF means also a 
larger Qa- As a consequence the B production in MF is 
larger and, because these B's have to be able to leave the 
pipe via desorption, the probabilities (Bi) and (B5) are 
larger in MF. The probabilities (Ai) and (A5) are there- 
fore smaller, which means that the MF curves and the 
simulation curves in figure cross each other, as can actu- 
ally be seen. The behavior of the system at high loading 
and at low loading is about the same, except that A is 
smaller at high loading. 

One might expect that the larger the number of reactive 
sites the more B's will be produced in the pipe. From the 
simulations we see that the amount of B's produced per 
unit time by all reactive sites goes to a limit value when 
the number of reactive sites is increased. In figure |, the 
marked line represents the B production as a function 
on the length of the pipe and the dashed line the B pro- 
duction according to MF. For short pipe lengths, the B 
production from both MF and simulations increase lin- 
early with S, while for higher lengths it converges to a 
limiting value. The limiting value is higher for MF. This 
could be seen also from the Table ||. According to MF 
there are more B produced in the pipe. 



00 we have 



From the simulations (see figure ^ we see that for high 
adsorption rates, Bprod converges to a point and the 
corresponding value is equal to the analytical value for 
the case adsorption is infinitely fast. The reason for 
this is that all the sites are occupied, diffusion is com- 
pletely suppressed, and only the marginal sites play a 
role. The expression above can be seen as a factor 2 
for the two marginal sites, the probability that an A at 
the marginal sites is converted to a B before it desorbs 
Wrx/VFrx + Wdcs, and the rate constant for desorption 

The accuracy of the simulation results for Qa and 
Bprod can be derived by looking at the total loading Q in 
Table ||. For the total loading Q, the simulation results 
can be compared with the values of the exact expression 
(12). We remark that the largest deviation from the ex- 
act analytical results is 0.04, so the relative errors are 
around 0.02%. 

The differences between MF and the simulations be- 
comes especially clear in the limit Wdig 00. Because 
this limit makes the system homogeneous in MF we get 



Qi 



ads 



w, 



ads 



o 



ads 



Wads + Wdcs Wrx + Wdc 



(38) 



(39) 



The first factor in these expressions is the probability 
that a site is occupied. The second factor indicates if the 
particle is converted to a B or not before it desorbs. The 
simulations show that the system should not be homoge- 
neous at all (see figure ^). The B production increases 
linearly with S only for the case of infinitely-fast diffu- 
sion, otherwise it converges to a limiting value. 



C. Only some of the sites reactive 

We consider now the situation that not all the sites 
are reactive, and that these reactive sites can be either 
uniformly distributed inside the pipe or distributed in 
compact blocks. We will show that the number of reac- 
tive sites doesn't change qualitatively the properties of 
the system. Qa, Qb and number of B produced for a 
variable number of reactive sites are compared with the 
previous results. 



1. Mean Field 

From the Master Equation it is easy to show that the 
total loading is again just the same as in the case when 




FIG. 4: The site occupancy with A ((A„)) and B {{B„)) as a function on site number for cases a, b, c, d, e when S = Nreac = 30. 
The continuous Une and the corresponding symmetric Une represent MF results for (A„) and (B„) respectively. The dashed 
lines represent DMC results for {A„) and {B„). (A„) is decreasing towards the middle of the pipe while (B„) is increasing. 



0.02 



I ' ' ' ' ' ' ' ' ' 1 

5 10 15 20 25 30 35 40 45 50 

S 



FIG. 5: B production as a function on the length of the pipe for Wads=0.2, Wdes=0.8, Wdiff=2, and Wrx=0.1. The marked 
line represents the DMC results and the dashed line represents the MF results. 




FIG. 6: Analytical and simulation results for site occupancy of a system when parameters are: S — Nreac = 30, Wads=0.8, 
W^des=0.2, Wdiff=100 and Wrx=2. The continuous line and the corresponding symmetric line represent the simulation profiles 
for site occupancy with A and B particles. The bottom and the upper straight lines represent the analytical results for 
occupancy with A and with B particles respectively. 



all the sites are reactive. We introduce an extra coef- 
ficient A„ in the MF equations to the reaction term. 
A„ = 1 if n is a reactive site and A„ = if it is not 
a reactive site. The steady-state equations are identical 
to equations (25), except that Wrx should be replaced 
by WrxA„. The resulting set of equations is linear again 
and it should be possible to solve them numerically. In 
fact only the probabilities for the marginal and reactive 
sites have to be solved numerically. For the other sites 
the probabilities can be obtained by simple linear inter- 
polation. That this is correct can be seen because those 
sites only have the diffusion term. We can also remove 
the probabilities for the B's, because we have from the 
model without conversion that 



(A„) + (B„) = 1 - (*„) = 



ads 



(40) 



The resulting equations for the reactive sites have the 
same form as equation (25) for the non-marginal sites. 
We expect therefore that we get an exponential decrease 
of (A„) on the reactive sites when we move from the 
marginal sites to the center of the pipe, and a linear de- 
pendence on n between the unreactive sites. 



2. 



Simulation results 



The number of reactive sites is considered to vary from 
1 to 50% and the reactive sites are distributed either in 
blocks situated near the marginal sites, in the middle of 
the pipe, or homogeneously distributed in the pipe. We 
will first compare the MF results with the MC simula- 
tion results for different sets of parameters and then we 
look at the dependence of B production and total load- 
ing Qa on the number and position of reactive sites. For 
the comparison between MF and MC results we consider 
the system size S = 30 and the number of reactive sites 
Nreac = 10. The scts of parameters used for the specific 
situations to be studied are the same as the sets used in 
the case with all the sites reactive in the previous section. 

We can see from the Tables [II and [V that when the 



reactive sites are homogeneously distributed or situated 
as a block in the middle of the pipe, there are significant 
differences between MF results and MC results. When 
the reactive sites form blocks near the marginal sites, the 
results are almost the same as when all sites are reactive: 
the MC and the MF results differ if we have fast diffusion 
and slow reaction. The sites in the center of the pipe are 
not relevant when the sites at the ends of the pipe are 
reactive. When the reactive sites are situated only in the 
middle of the pipe, we have deviations for all the sets of 
parameters. They are very prominent for the case when 
we have high loading, fast diffusion and fast reaction. 
MF strongly underestimates A's for all non-reactive sites, 
but we have also important deviations for high loading in 
the cases with fast diffusion-slow reaction, slow diffusion- 
fast reaction, and slow diffusion-slow reaction. This is 
happening because for high loading, the end sites will 



always be occupied by a particle A and the B's will not 
be able to get out of the pipe. In MF particles effectively 
can pass each other, so B particles are then able to get out 
of the pipe. Even for the case of low loading we still have 
deviations from MF for fast diffusion and fast reaction. 
In this case MF overestimates A's for nonreactive sites. 
For fast diffusion and slow reaction, MF underestimates 
A's for nonreactive sites and for slow diffusion and slow 
reaction MF overestimates B's for the reactive sites in 
the middle. 

Figures ^ and || show how the probabilities (A„) and 
(B„) vary in the pipe. The situations for reactive sites 
forming blocks at the ends of the pipe are not shown 
as they are almost the same as when all the sites are 
reactive (see figure 4). When the reactive sites are ho- 
mogeneously distributed the plots look also very similar 
to the ones with all sites reactive, except that the char- 
acteristic length A is larger. (A„) and (B„) look very 
different when the reactive sites form a block in the mid- 
dle of the pipe. The MF results show, as predicted, a 
linear behavior at the nonreactive sites. The MC results 
show, however, a nonlinear behavior in the form of S-like 
curves. At the reactive sites the behavior is similar to 
the situation with all sites reactive with the MC results 
showing a more rapid approach to the value at the mid- 
dle of the pipe than MF, i.e., smaller A. The values at 
the marginal sites can differ between MC and MF quite 
a lot. This reflects the difference in Bprod mentioned be- 
fore: a different Bprod must be accompanied by a different 
B desorption at steady-state. As we have already seen 
from the case when all the sites were reactive, Bprod very 
rapidly approaches the limiting value when the pipe is 
made longer (see figure 5). Similarly when we start with 
few reactive sites and, instead of increasing the length of 
the pipe, we increase the number of reactive sites. The 
loading Qb is already almost the same as the value with 
all sites reactive when only about 10% of all sites are re- 
active provided there are reactive sites at or very near 
the marginal sites. If the reactive sites are moved away 
from the ends of the pipe, then the loading Qb and the 
B production decreases. 



IV. 



SUMMARY 



We have used analytical and simulation techniques to 
study the reactivity in Single-File Systems. 

The MF results show that MF models Single-File be- 
havior by changing the diffusion rate constant, but it 
effectively does allow passing of particles. 

When all the sites are reactive, the simulation and MF 
results are very similar for all the parameters, except 
for the case when we have low reaction rates and fast 
diffusion. In these cases MF overestimates the amount of 
A's in the pipe. The amount of B produced per unit time 
by all reactive sites goes to a limit value when the number 
of reactive sites is increased. For high adsorption rates, 
Bprod converges to a point and the corresponding value 





Marginal 




Middle 


Homogeneous 


set 


MF 


Sim 


MF 


Sim 


MF 


Sim 


a) 


0.0512 


0.0500 


0.0731 


0.0771 


0.0208 


0.0469 


b) 


0.0153 


0.0152 


0.0712 


0.0788 


0.0138 


0.0206 


c) 


0.0590 


0.0672 


0.0901 


0.0881 


0.0719 


0.0594 


d) 


0.0041 


0.0041 


0.0667 


0.0730 


0.0120 


0.0123 


c) 


0.0067 


0.0006 


0.0447 


0.0583 


0.0126 


0.0121 


f) 


0.0896 


0.0752 


0.3008 


0.3473 


0.1121 


0.1056 


g) 


0.0376 


0.0369 


0.2850 


0.3383 


0.0585 


0.0605 


h) 


0.0871 


0.0579 


0.2844 


0.3250 


0.1227 


0.0867 


i) 


0.0048 


0.0048 


0.2606 


0.3137 


0.0319 


0.0413 


j) 


0.0056 


0.0053 


0.1556 


0.2826 


0.0175 


0.0289 



TABLE III: Simulation and MF results for Qa for all the sets of parameters in the cases of homogeneous distribution of the 
reactive sites, blocks of reactive sites in the middle of the pipe and near the marginal sites (5* = 30, Nreac = 10) 







Marginal 




Middle 


Homogeneous 


set 


MF 


Sim 


MF 


Sim 


MF 


Sim 


a) 


0.0099 


0.0099 


0.0001 


0.0000 


0.0011 


0.0045 


b) 


0.0449 


0.0477 


0.0001 


0.0008 


0.0153 


0.0117 


c) 


0.1156 


0.1021 


0.0393 


0.0216 


0.0728 


0.0521 


d) 


0.2449 


0.2492 


0.0283 


0.0161 


0.1357 


0.0891 


e) 


0.2767 


0.2763 


0.1482 


0.0646 


0.2410 


0.1899 


f) 


0.0239 


0.0235 


0.0014 


0.0006 


0.0086 


0.0076 


g) 


0.1129 


0.1160 


0.0015 


0.0000 


0.0139 


0.1171 


h) 


0.1796 


0.1421 


0.0470 


0.0059 


0.1212 


0.0661 


i) 


0.2931 


0.2041 


0.0288 


0.0069 


0.1526 


0.0897 


j) 


0.30i;i 


{):2<.m 


().i5."32 


0.0il3 


0.2713 


0.17:!iJ 



TABLE IV: Simulation and MF results for Bprod for all the sets of parameters in the cases of homogeneous distribution of the 
reactive sites, blocks of reactive sites in the middle of the pipe and near the marginal sites (5* = 30, Nreac = 10) 

is equal to the analytical value for the case adsorption is 
infinitely fast. The sites in the middle of the pipe have 
no effect on the B production. The differences between 
MF and the simulations becomes especially clear in the 
limit Wdiff oo. 

When only some of the sites are reactive, there are sig- 
nificant differences between MF and MC results when the 
reactive sites arc homogcmcously distributed or situated 
as a block in the middle of the pipe. When the reactive 
sites form blocks near the marginal sites, the results are 
almost the same as when all sites are reactive: The MC 
and the MF results differ only when we have fast diffu- 
sion and slow reaction. The sites in the center of the pipe 
arc not relevant when the sites at the ends of the pipe 
are reactive. When the reactive sites are situated in the 
middle of the pipe, we have deviations for all the sets of 
parameters. They are very prominent for the case when 
we have high loading, fast diffusion and fast reaction. MF 
strongly underestimates A's for all non-reactive sites, but 
we have also important deviations for high loading in the 
cases with fast diffusion-slow reaction, slow diffusion-fast 
reaction, slow diffusion-slow reaction. The MF results 
show a linear behavior at the nonreactive sites. The MC 
results show, however, a nonlinear behavior in the form 
of S-like curves. The loading Qb is already almost the 
same as the value with all sites reactive when only about 



10% of all sites are reactive provided there are reactive 
sites at or very near the marginal sites. If the reactive 
sites are moved away from the ends of the pipe, then the 
loading Qb and the B production decreases. 
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VI. APPENDIX 

A. Probability to find the system in a certain 
conflguration. Loadings and fluctuations. 

We show the existence of a function depending only 
on the number of particles such that 

P„ = q{n{a)) (41) 

is the steady-state solution of the Master Equation (2) for 
a system without conversion, where n{a) is the number 
of particles in configuration a. The second part of the 
proof consists of showing the uniqueness of the solution. 
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FIG. 7: The site occupancy for the cases a, b, c, d, e -homogeneous distribution. The continuous line and the corresponding 
symmetric line represent the MF results. The other dashed lines represent the DMC results. 
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FIG. 8: The site occupancy for the cases f, g, h, i, j - middle sites reactive, Nreac = 10. The continuous line and the 
corresponding symmetric line represent the MF results. The other dashed lines represent the DMC results. 



Substitution of Pa = q{n(a)) in equation (9) shows 
that the last term in the Master Equation vanishes, be- 

and n{a) = The other terms 



cause A 



(diff) 

a [3 



(diff-) 



can also be simplified by using how the number of parti- 
cles changes upon adsorption and desorption. 



at 



5(n(a)-l)^Ai^;^)-,(n(a))^A^ 



(ads) 
13a 



fdes) 



+ Wdes 5(n(a) + 1) ^ Ay - q{n{a)) ^ Ag 

(42) 

A further simplification is possible if we realize that des- 



orption reverses the effect of an adsorption and vice versa. 
This means A^')^""^ = This leads to 

^ = [q{n{a) - l)W^. - <z(n(a)) W^eJ E^t^ 

- [g(n(a))W^ads - q{n{a) + l)T4-des] E^t'^' 

(43) 

We denote by N, the number of particles in a certain 
configuration, A'' = n{a). We see that we get a steady- 
state solution for 



dt 



0, 



(44) 



provided by 



q{N + l) ^ Wad. 
q{N) Wdc. 



(45) 



for TV = 0, 1, 2, ... , S-1. (Note that the case N ^ S in 
the Master Equation presents no problems, because the 
summation over /3 yields zero.) 

The second step consists of showing that this solution 
is the only one. This part for instance can be found in 
Chapter 5 of Van Kampen. [pl| 



B. Derivation of function q(N) 

Expression (45) leads to 
q{N)=C 



Wade 



Wdc 



(46) 



where C is some normalization constant. We can com- 
pute it from 

5 



N=0 



Wads 



Wdc 



N 



N=0 



c 



q{N) 



Wdcs + Wads 



Wde 



(47) 



The combinatorial factor after the third equal sign de- 
rives from the number of configurations with N particles. 
The last step uses 



The expression for q{N) now becomes 
qiN) 



Wdes 


S 


[Wads] 




Wdcs + Wads. 




[ Wdcs J 





(48) 



(49) 



Note that this expression does not depend on Wdig: i.e., 
diffusion has no effect at all on steady-state properties. 

The probability p{N) that there are N particles in the 
system is given by 



p{N) 



q{N) 



Wde 



Wdes + W, 



ads 



w, 



ads 



Wdc 



(50) 



This follows from (49). With this formula we can com- 
pute all statistical properties of the number of particles. 
The average number of particles is 



E ^p(^^ = 



N=0 



Wd, 



Wdes + Wads 
Wads 



s s 

E 



N 



ads 



Wde 



N 



Wdes + Wads 



s. 



The loading of the pipe, defined as the average number 
of particles per site, is 



Wads 



S Wads + Wde 

The average squared number of particles is 

n 5 S 



(52) 



N=0 



Wde 



Wdes + Wads 



E 



Wads(Wdes + ^Wads) 
(Wdes + Wads)2 



W, 



n N 



ads 



Wde 



s. 



(53) 



The variance, i.e., the square of the fiuctuation in the 
number of particles, is then 



Y,N'p{N) 



N=0 



E ^p^^) 



.N=0 



Wads Wde 



(Wdes + Wads 



r5. 



(54) 



by 



C. Derivation of the one-site and two-sites 
occupancy for the model without conversion 

The probability that site n is occupied by A is given 



(A„)==^P„A("), 

a 

= EE^"^^"^ 

Af QGAf 

= E9(^)E^^"^ 



N 
S 



E9(^) 



Af=l 



Wde 



ctGAf 

5-1 
N -l^ 

s s 



(55) 



E 

Af=l 





■Wads" 




[ Wdes J 



Af 



w. 



ads 



Wde 



ads 



Wde 



S-1 



(51) 



Wdes + Wads 

Wdes 
Wdes + Wads 
^ Wads 

Wads + Wdes ' 

where aL""* is 1 if site n in configuration a is occupied by 
an A particle, and it is otherwise. The combinatorial 
factor denotes the number of ways the particles except 
the one at site n can be distributed over the remaining 
sites. Knowing the one-site occupancy we can derive the 
two-site occupancy 



D. Continuum limit 



(A„A„+i) = 5^P„Ai")Ai"+i), 

a 



= E ( 



Wde 



5-2' 
N-2j ' 

s s 



W^des 



E 

N=2 



Wdes + Wads, 
Wads + Wdes 



^ 'Wads' 



w 



des 



1 + 



w 



ads 



dcs 



AT 



S-2 



(56) 



The rate equation for the A's is 

- (*„_iA„)] - Wrx(A„). 



= Wdifi[(A„_i*„) + (*„A„+i) - (A„*„+i) ^g^^ 



The MF approximation of this equation is 
d(A„) 



dt 



Wdiff[((A„_i) + (A„+i))(*„) 

(A„)((*„_i) + (*„+!))] (58) 
Wx(A„). 



If we take the continmini Hmit and denote by a = a(x, t), 
b = b{x,t) and v = v{x,t) the probabiUty distribution of 
A's, B's and vacancies respectively, and if we use Taylor 
series for the diffusion term, the equation becomes 



= Wdiff [a((n - l)d, t) + a{{n + l)d, t)]v{nd, t) 



WiiQa{nd,t)[v{{n — l)d,t) + v{{n + l)d,t)]a{nd,t) — W^a 
^ Wdiff 



,9a 1 ,^ d'^a ,da 1 9^a"' 

ox 2 a.T^ ox 2 ax^ 



- Wdifi 
= Wdiff 
= Wdiffd^ 
= Wdiffd' 



dx 2 2 dx'^ 



2a + d 



dx^ 



V - Wdifffl 



2v + d 



dx"^ dx"^ dx^ 



dx^ 



a — WrxCt 
- Wrxfl 



- Wrxfl, 



r 



(59) 



where d is the distance between sites. A similar relation 
can be derived for b{x, t). With D = Wdifirf^ we can write 

I 



f da/dt\ _ jj f ^ — b a \ f d'^a/dx 



\db/dt ) 



r 



(60) 



E. MF derivation of the total loading Qa in case so, the expression for Qa is 
with conversion 

Qa = |E[^i"+^i^^'""]' 

n=l 
S 



The total loading with A's, Qa, is written as 

s 

S 



(61) 



Xi 



S+l-n 



n=l 

S 



71=1 



(62) 
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Steady-State Properties of Single-File Systems with Conversion 



S.V. Nedea,* A.P.J. Jansen,t J.J. Lukkien,t and P.A.J. Hilbers^ 
(Dated: January 31, 2002) 

Wo have used Monto-Carlo methods and analytical techniques to investigate the influence of the 
characteristic parameters, such as pipe length, diffusion, adsorption, desorption and reaction rate 
constants on the steady-state properties of Single-File Systems with a reaction. We looked at cases 
when all the sites are reactive and when only some of them are reactive. Comparisons between 
Mean-Field predictions and Monte Carlo simulations for the occupancy profiles and reactivity are 
made. Substantial differences between Mean-Field and the simulations are found when rates of 
diffusion are high. Mean-Field results only include Single-File behavior by changing the diffusion 
rate constant, but it effectively allows passing of particles. Reactivity converges to a limit value 
if more reactive sites are added: sites in the middle of the system have little or no effect on the 
kinetics. Occupancy profiles show approximately exponential behavior from the ends to the middle 
of the system. 

PACS numbers: 02.70Uu, 02.60.-x, 05.50.-|-q, 07.05.Tp 



I. INTRODUCTION 

Molecular sieves are crystalline materials with open 
framework structures. Of the almost two billion poimds 
of molecular sieves produced in the last decade, 1.4 billion 
pounds were used in detergents, 160 millions pounds as 
catalysts and about 70 millions pounds as adsorbents or 
desiccants. [1] 

Zeolites represent a large fraction of known molecular 
sieves. These are all aluminosilicates with well-defined 
pore structures. In these crystalline materials, the metal 
atoms (classically, silicon or aluminum) are surrounded 
by four oxygen anions to form an approximate tetrahe- 
dron. These tetrahedra then stack in regular arrays such 
that channels and cages are formed. The possible ways 
for the stacking to occur is virtually unlimited, and hun- 
dreds of unique structures are known. [2] 

The channels (or pores) of zeolites generally have cross 
section somewhat larger than a benzene molecule. Some 
zeolites have one-dimensional channels parallel to one an- 
other and no connecting cages large enough for guest 
molecules to cross from one channel to the next. The 
one-dimensional nature leads to extraordinary effects on 
the kinetic properties of these materials. Molecules move 
in a concerted fashion, as they are unable to pass each 
other in the channels. These structures are modeled 
by one-dimensional systems called Single-File Systems 
where particles are not able to pass each other. A par- 
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tide can only move to an adjacent site if that site is not 
occupied. 

This process of Single-File diffusion has different char- 
acteristics than ordinary diffusion which affects the na- 
ture of both transport and conversion by chemical re- 
actions. For Single-File diffusion, the mean-square dis- 
placement of a particular particle is proportional to the 
square-root of time 

(r^) = 2Fti 

where F is the Single-File mobility. [3] This is in contrast 
to normal diffusion, where mean-square displacement is 
directly proportional to time. A variety of approaches 
have been used to describe the movement of the particles 
in Single-File Systems, most of them concentrated on the 
role of the Single-File diffusion process. 

Molecular Dynamic(MD) studies of diffusion in zeo- 
lites have become increasingh^ popular with the advent of 
powerful computers and improved algorithms. In a MD 
simulation the movement is calculated by computing all 
forces exerted upon the individual particles. MD results 
have been found to match experimental observations of 
Single-File diffusion for systems with one type of molecule 
without conversion and with very short pores. [4-7] Be- 
cause a molecule can move to the right or to the left 
neighboring site only if this site is free, MD simulations 
under heavy load circumstances require a high computa- 
tional effort for particles that hardly move. However, the 
level of detail provided by MD simulations is not always 
necessary. 

Thus, deterministic models are used also but they are 
mainly focused on dynamic and steady-state information 
of short pore systems. [8-10] Several researchers [11- 
13] used a stochastic approach, i.e.. Dynamic Monte 
Carlo(DMC), to determine the properties of Single-File 
Systems. In DMC reactions can be included. The rates 
of the reactions determine the probability with which dif- 
ferent configurations are generated and how fast (at what 
moment in time) new configurations are generated. The 
most severe limitation of the DMC method arises when 



the reaction types in a model can be partitioned into 2 
classes with vastly different reaction rates. In this case, 
extremly large amounts of computer time are required 
to simulate a reasonable number of chemical reactions. 
However, in general the system can be simulated for much 
longer times than with MD. 

All the previous references put the emphasis on the 
transport properties of adsorbed molecules as the impor- 
tant factor in separation and reaction processes that take 
place within zeolites and other shape-selective microp- 
orous catalysts. Rodenbeck and Karger [9] solved numer- 
ically the principal dependence of steady-state properties 
such as concentration profiles and the residence time dis- 
tribution of the particles, on the system parameters for 
sufficiently short pores. In multiple papers, Auerbach et 
al. [14, 15] used Dynamic Monte Carlo to show different 
predictions about Single-File transport and direct mea- 
surements of intercage hopping ion strongly adsorbing 
quest-zeolite systems. Saravanan and Auerbach [16, 17] 
studied a lattice model of self-diffusion in nanopores, to 
explore the influence of loading, temperature and adsor- 
bate coupling on benzene self-diffusion in Na-X and Na-Y 
zeolites. They applied Mean-Field (MF) approximation 
for a wide set of parameters, and derived an analyti- 
cal diffusion theory to calculate diffusion coefficients for 
various loadings at fixed temperature, denoted as "diffu- 
sion isotherms" . They found that diffusion isotherms can 
be segregated into subcritical and supercritical regimes, 
depending upon the system temperature relative to the 
critical temperature of the confined fiuid. Supercritical 
systems exhibit three characteristic loading dependen- 
cies of diffusion depending on the degree of degeneracy 
of the lattice while the subcritical diffusion systems are 
dominated by cluster formation. Coppens and Bell [18- 
20] studied the influence of occupancy and pore network 
topology on tracer and transport diffusion in zeolites. 
They found that diffusion in zeolites strongly depends 
on the pore network topology and on the types and frac- 
tions of the different adsorption sites. MF calculations 
can quickly estimate the diffusivity, although large devi- 
ations from the DMC values occur when long-time corre- 
lations are present at higher occupancies, when the site 
distribution is strongly heterogeneous and the connectiv- 
ity of the network low. 

Few researchers included also reactivity in Single-File 
Systems. Tsikoyannis and Wei [8] considered a reac- 
tive one-dimensional system with all the sites reactive 
in order to get more information about the reactivity 
and selectivity in one-dimensional systems. They used a 
Markov pure jump processes approach to model zeolitic 
diffusion and reaction as a sequence of elementary jump 
events taking place in a flnite periodic lattice. Monte 
Carlo and approximate analytical solutions to the de- 
rived Master Equation were developed to examine the 
effect of intracrystalline occupancy on the macroscopic 
diffusional behavior of the system. One conclusion was 
that better results using analytical approach can be ob- 
tained compared to DMC simulation results by including 



more correlations between neighboring sites in regions of 
the systems with high occupancy gradients and less cor- 
relations in regions with low and no occupancy gradients. 
Starting from Wei [8] results about correlations in Single- 
File Systems, Okino and Snurr [10] used a deterministic 
model where each site was assumed to have equal activ- 
ity towards reaction. Doublet approximation was found 
to overpredict the occupancy of the sites and the increas- 
ing mobility raised the concentration of reactants in the 
pore. 

Using DMC simulations we have observed that even for 
inflnitely fast diffusion, we still have Single-File effects in 
the system. Instead of focusing on diffusion at different 
occupancies of the system, we therefore concentrate in 
this paper on the reactivity of the system, studying the 
reactivity of the system for different sets of kinetic pa- 
rameters, the length of the pipe and the distribution of 
the reactive sites. We analyse the situations when MF 
gives good results and when MF results deviate strongly 
from the DMC simulations. We investigate the effect of 
the various model assumptions made about diffusion, ad- 
sorption/desorption, and reaction on the overall behav- 
ior of the system. We look at the total loading, loading 
with different components, generation of reaction prod- 
ucts and occupancies of individual sites as a function of 
the various parameters of a Single-File System. 

In section II we specify our mathematical model for 
diffusion and reaction in zeolites together with the theo- 
retical background for the analytical and simulation re- 
sults. In section III A we present the various results for 
the simplifled model without conversion. In section IIIB 
we use MF theory to solve the Master Equation govern- 
ing the system behavior for the case when all the sites 
have the same activity towards conversion. Similarly the 
results obtained using DMC simulations are presented in 
section IIIB 2 and are compared with MF results. We 
pay special attention to the infinitely fast diffusion case 
and to the infiuence of the length of the pipe on the over- 
all behavior of the system. In section IIIC we analyze 
again the MF and simulation results but for the case 
when only some of the sites are reactive. The infiuence 
of the position and number of reactive sites on the reac- 
tivity and site occupancy of the system is outHned. The 
last section summarizes our main conclusions. 



II. THEORY 

In this section we will give the theoretical background 
for our analytical and simulation results. First we wfll 
specify our model and then we will show that the defined 
system obeys a Master Equation. [21] We will simulate 
the system governed by this Master Equation using DMC 
simulations. The rate equations used for the derivation 
of the analytical results are outlined. 



A. The Model 



B. Master Equation 



Because we are interested in reaction of molecules in 
Single-File Systems, we call the system we are modelling, 
Single-File System with conversion. We model a Single- 
File System by a onc-dimensional array of sites, each 
possibly occupied by a single adsorbate. The sites are 
numbered 1, 2, . . . , S. An adsorbate can only move if 
an adjacent site is unoccupied. The sites could be reac- 
tive or unrcactivc and wc note with N,.eac the number of 
reactive sites. A reactive site is the only place where a 
reaction may take place. 

Wc consider two types of adsorbatcs, A and B, in 
our model and we denote with X the site occupation 
of a site, X=(*, A, B), which stands for an empty 
site, a site occupied by A, or a site occupied by a B, 
respectively. The sites at the ends of the system are 
labeled with m, and the reactive sites are labeled with 
r (see figure 1). We restrict ourselves to the following 
mono and bi-molccular transitions. 

a) Adsorption and desorption 

Adsorption and desorption take place only at the 
two marginal sites i.e., the left and rightmost sites at 
the ends of the system. 



A(gas) -I- — > Am 
Am — > A(gas) + *m 
Bm — > B(gas) -I- *m, 

where m denotes a marginal site. Note that there is no 
B adsorption. B's are formed only by a reaction. 



b) Diffusion 

In the pipe, particles are allowed to diffuse via hopping 
to vacant nearest neighbor sites. 



A„ -I- *n+l < > *n + 

Bn + *n+l ■* * *n + B„+j , 

where the subscripts are site indices: n=l, 2, . . . , S-1. 
c) Reaction 

An A can transform into a B at a reactive site. 



Reaction kinetics is described by a stochastic process. 
Every reaction has a microscopic rate constant associated 
with it that is the probability per unit time that the 
reaction occurs. Stochastic models of physical systems 
can be modelled by a Master Equation. [21] 

By a, (5, we will indicate a particular configuration of 
the system i.e., a particular way to distribute adsorbates 
over all the sites. Fait) will indicate the probability of 
finding the system in configuration a at time t and Wai3 is 
the rate constant of the reaction changing configuration 
to configuration a. 

The probability of the system being in configuration 
a at time t + dt can be expressed as the sum of two 
terms. The first term is the probability to find the sys- 
tem already in configuration a at time t multiplied by 
the probability to stay in this configuration during dt. 
The second term is the probability to find the system in 
some other configuration /3 at time t multiplied by the 
probability to go from P to a during dt. 

Pa{t + dt) = il-dtJ2W0c.)Pait)+dtY,WapPpit) (1) 

P P 

By taking the Hmit dt ^ this equation reduces to a 
Master Equation: 



dPajt) 

dt 



E[W^a/3^/3W-W^/3aJ'aW]. (2) 

P 



Analytical results can be derived as follow. The value 
of a property X is a weighted average over the values 
which is the value of X in configuration a: 



{X)=J2PaXa- (3) 

a 

From this follows the rate equation 

d{X) _^dP^^ 
dt ~ ^ dt " 

a 

= ^[^a/3f/3 - WpaPa]Xa (4) 
aP 

= J2Wc.pPp{Xc,-Xp). 

aP 



The initial state of the system is all that all sites are 
empty (no particles in the pipe). In this paper we will 
only look at steady-state properties and not to the time 
dependence of the system properties starting with no par- 
ticles. 



C. Dynamic Monte Carlo 

Because it might be not always possible to solve the 
Master Equation analytically, DMC methods allow us to 
simulate the system governed by the Master Equation 
over time. We simplify the notation of the Master Equa- 
tion by defining a matrix W containing the rate constants 




I adsorption 




FIG. 1: Picture of a Single-File System with two types of adsorbates, A(lighter colored) and B(darker colored). The marginal 
sites are labeled with m, and the reactive sites(lighter colored) with r. Adsoption of A and desorption of A and B can take 
place only at the two marginal sites. An A can transform into a B only on r labeled sites. 



ap, and a diagonal matrix R by Rap = ^^iP^ 



, W^p, if 

a = P, and otherwise. If we put the probahiHties of 
the configurations Pa in a vector P, we can write the 
Master Equation as 



dP 

dt 



-(R- W)P. 



(5) 



where R and W are assumed to be time independent. 
We also introduce a new matrix Q, Q(i) = exp[— Ri]. 

This matrix is time dependent by definition and we 
can rewrite the Master Equation in the integral form 



P{t) = Q{t)P{0)+ f dt'Q{t~t')WP{t'). (6) 
Jo 

By substitution we get of the right-hand-side for P{t') 



p{t) = m) 

+ [ dt'Cl{t-t')WCl{t') 
Jo 



ft rt' 

+ 

lo Jo 
+ ...]P(0). 



/ dt' [ dt"Q{t - t')WQ{t' - t")WQ{t" 
Jo Jo 



(7) 



Suppose at i = the system is in configuration a 
with probability Pa(0). The probability that, at time 
t, the system is still in configuration a is given by 
Qaa{t)Pa{0) = Raat) Pa{0) ■ This shows that the 

first term represents the contribution to the probabilities 
when no reaction takes place up to time t. The ma- 
trix W determines how the probabilities change when 
a reaction takes place. The second term represents the 
contribution to the probabilities when no reaction takes 
place between times and t' , some reaction takes place at 
time t' , and then no reaction takes place between t' and t. 
The subsequent terms represent contributions when two, 
three, four, etc. reactions take place. The idea of the 
DMC method is not to compute probabilities Pa{t) ex- 
plicitly, but to start with some particular configuration. 



representative for the initial state of the experiment one 
wants to simulate, and then generate a sequence of other 
configurations with the correct probability. The method 
generates a time t' when the first reaction occurs accord- 
ing to the probability distribution 1 — exp[—Raat\- At 
time t' a reaction takes place such that a new configura- 
tion a' is generated by picking it out of all possible new 
configurations l3 with a probability proportional to Wa'a- 
At this point we can proceed by repeating the previous 
steps, drawing again a time for a new reaction and a new 
configuration. 



III. RESULTS AND DISCUSSION 
A. No conversion 

We mention in this section various results for the sys- 
tem without conversion. These results can be derived 
analytically. The derivations are not difficult, so for com- 
pleteness we give them in the appendix. We will use the 
results when we deal with the system with conversion. 

In a Single-File System without conversion, the rele- 
vant processes to describe are adsorption, desorption and 
diffusion. So, Wap is given by 



Wap = Wad.Ai7' + T^de.Ai7' + T^difFA^ , 



(8) 



where A^'^^ equals 1 if a reaction of type "rx" can trans- 
form the system from /] to a, and equals otherwise. 
W^ads, W^des, Wdiff are the rate constants of adsorption, 
desorption and diffusion respectively. 

If we substitute expression (8) into the Master Equa- 
tion (2), we get 



dPg 

dt 



Wdiff Yl 



(9) 



-I- 



(difF) p 

'aP ^13 



lifF) p 1 



Using this expression we can show that when the sys- 
tem is in steady state then the probabihty of Ending the 
system in a certain conGguration depends only on the 
number of particles in the system. 



(10) 



where is the number of particles in configuration a. 
The expression for q{N) is: 



qiN) 





S 


W'ads 


N 


W^des + W^ads. 




[W^desj 





(11) 



Note that diffusion has here no effect on steady-state 
properties. 

The loading of the pipe, defined as the average number 
of particles per site, is then 



1 ^ 



Wads 



JV=0 



Wads + Wdes' 



(12) 



where p{N) is the probability that there are N particles 
in the system. Note again that diffusion doesn't influence 
the steady-state loading. 

The standard deviation, i.e., the fluctuation in the 
number of particles is then: 



(13) 




(Wdes + Wads) 



To determine how the parameters of the system influ- 
ence the kinetics of the system, we are interested in the 
correlation in the occupancy between neighboring sites. 
We look at one site occupancy and at two sites occupan- 
cies. We denote by (A„) the probability that an A is at 
site n and with (A„A„_|_i) the probabflity to have an A 
at site n and one at site n + 1. 

One and two-site probabilities can be derived from the 
fact that all conflgurations with the same number of par- 
ticles have equal probabflity and the expressions for q{N). 
We find 



and 



(A„) 



(A„A„+i) 



Wads 



Wads + Wdes' 



Wads 



Wads + Wdes 



(14) 



(15) 



Note that this probability does not depend on the site, 
all sites have equal probability to be occupied and that 
there is no correlation between the occupation of neigh- 
boring sites. Again diffusion doesn't influence these prop- 
erties. Note also that these expressions are the same as 
for a model in which particles are aflowed to pass each 
other. 



B. All sites reactive 

We look flrst at the situation with all sites reac- 
tive: i.e., conversion of an A into a B particle can 
take place at any site including the marginal sites. For 
simplicity we consider WciesA=WdesB=Wdes, and also 
WdifFA=WdifrB=Wdiff. We wifl be looking at the total 
loading (Q), the total loading of A's (Qa), the total load- 
ing of B's (Qb), the number of B's produced per unit time 
(Dprod), and how the distribution of A's and B's varies 
from site to site ((A„) and (B„)). 

Note that the total loading of the pipe for the model 
with conversion is the same as for the model without 
conversion 



Q 



Wads 



Wads + Wde 



(16) 



The loadings and the production of B's can easily be 
derived from the probabilities {A„) and (B„) so we flrst 
focus on them. For a non-marginal site we can write 



d{A„} 
dt 



D (A,diff) , D(rx) 



(17) 



where i?„^*''^'*^ is the rate of diffusion of A from and 
to site n, and RrS^^^ is the rate of conversion of A to 
B on site n. The conversion takes place at one site and 
is therefore easier to handle than the diffusion. Using 
equation (4) we have 



(rx) 



Wrx Aa^p 



(rx) 



-P/3(A„a — A„/3), 



(18) 



where A„„ = 1 if site n is occupied by an A in conflgu- 
ration a and A„a = if not. We have A„„ — A„^ 7^ if 
there is an A at site n in conflguration /3 (A„^ = l) that 
has reacted to a B leading to conflguration a (A„a = 0). 
This gives us 

= -W,^Y. = -WrAK), (19) 

where the prime restricts the summation to those /3's 
with Anp = 1. For the diffusion we simflarly get 

i?„(A>diff) = T^dis^A„/^'diff)p^(A„„ - A„^). (20) 

There are four ways in which A^^ — A„g 7^ and 
Aa^^*'^'"^' 7^ in /?: there is an A at site n that can move 
to site (n — 1), there is an A at site n that can move to 
(n + l), there is an A at site (n— 1) that can move to site 
n and there is an A at site (n+l) that can move to site n. 
In all cases we have Aa/g^^'*^'*' = 1. In the flrst two cases 
we have A„q — A„3 = — 1 and in the last two we have 
A„„ — A„^ = 1. The summation over j3 in the first case is 
restricted to configurations with an A at site n and a va- 
cant site (n — 1). This gives a term — WdifF(*n-i A„). The 



other cases give terms -VFdis (A„*„+i), Wdiff(A„_i*„), 
and Wdis{*n^n+i)- The rate equations then becomes 



dt 



Wdiff[-(A„*„_i) - (*„_iA„) + (A„_i*„) 



)] - W'rx(A„). 



(21) 



For (B„) we get similarly 



dt 



Wdiff[-(B„*„_i) - (*„_iB„) + (B„_i*„) 



+ (*„B„+i)] + PFrx(A„). 



(22) 



The marginal sites have also adsorption and desorption. 
They can be dealt with as the conversion. The rate equa- 
tions for A are 



rf(Ai) 
dt 

d{As) 
dt 



■■ Wdiff[-(Ai*2) + (*lA2)] + VFads(*l) 
Wdes(Ai)-W„(Ai), 

: Wdia[~{^S*S-l) + (*sAs-i)] + T^ads(*s) 



and the rate equations for B 



(23) 



rf(Bi) 

dt 
d{Bs) 

dt 



W^diff[-(Bl*2) + (*lB2)] - W^des(Bl) + W^rx(Ai) 
W^diff[-(As*S-l) + (*sAs-l)] - T^des(Bl) 



+ W^rx(As). 



(24) 



Note that these coupled sets of differential equations are 
exact. 



1. Mean Field results 

We will now look at the loadings Qa and Qb and the 
site occupation probabilities (A„) and (B„). We will 
first determine steady-state properties using the (MF) 
approximation: i.e, we put (A„*„_|_i) = (A„)(*„_|_i) etc. 
in the rate equations. This gives us 



Wads + Wdes 

WdisWdes 
Wads + Wdes 



+ 



Wads + Wde 
Wads Wdes 

Wads + Wdes ' 
WdiffWdes 

Wads + Wdes 
WdiffWd, 



[(A„+i) + (A„_i) - 2(A„)] - Wx(A„) 
[(B„+i) + (B„_i) - 2(B„)] + W,x(A„), 

[(A2) - (Ai)] - Wrx(Ai) - PFdes(Ai) 
[(B2) - (Bi)] + Wx(Ai) - Wdes(Bi), 



+ 



W'ads + W'des 
Wads Wdes 

Wads + Wdes ' 
WdiffH des 

Wads + Wdes 



^ [(As-i) - (As)] - W,x(As) - Wdes (As 



[(Bs-i) - (Bs)] + W,x(As) - Wdes(Bs 

(25) 



We have used here the probability for a site to be vacant 
that we have determined for the case without conversion. 

We note that these equations are identical to the MF 
equations of a system in which the particles can move 
independently with a rate constant for diffusion equal 
to WdiffWdes /(Wdes + Wads)- This means that the MF 
does not really model the non-passing that characterizes 
a Single-File System. 

The continuum limit of the MF equation is 



b 1-al \ d%/dx^ ' + 



(26) 

where a=a{x, t) is the probability distribution of A's (a 
similar definition holds for b), and D=Wdmd^ , with d the 
distance between neighboring sites (see appendix). These 
are the equations that arc normally used to describe dif- 
fusion in Single-File Systems. [19, 22-24] As this equation 
is derived from the MF equations, it has the same draw- 
back; i.e., the Single-File behavior is only incorporated 
by the reduction of the diffusion, but it does effectively 
allow for passing of particles. This shows up as so-called 
counter diffusion of A's and B's. [22-24] 

We see that equations (25) are linear and we can solve 
them at least numerically. We think however that it is 
worthwhile to use an analytical approach. We consider 
the ansatz 

(A„) cx a;" (27) 
in the steady-state equations (25) for (A„). This leads to 
- 2{l + a)x + l = Q, (28) 

with 

Wrx Wdes + Wads 



2Wdiff Wde 



(29) 
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FIG. 2: The characteristic length A as a function of xi. 



FIG. 3: Bprod per unit time at one marginal site as a function 

of Wads for S=Nreac=SO, W"des=0.8, WdilT = 2 and Wr,, = 0.4. 



The quadratic equation yields two solutions Xi and X2 
with X2 = x^^ . We have xi = X2 = ^ only when a = 0, 
i.e. when Wrx = 0. We will therefore assume a > and 
xi <1. Then 



xi = {l+a)- y/a{a + 2). 



(30) 



We can write then the solution (A„) = ai(a;i)" + 
02(3:1)'^"'"^"". The symmetry in the occupancy of the 
pipe (A„) = (As+i-n) yields ai = 02 = a. So, the gen- 
eral solution for the steady state has the form: 



(A„) =aK+xf+i-"). 



(31) 



The coefficient a is to be determined from the equa- 
tions for the marginal sites in the set of cquations(25). 
In the left side of the system n is small and {S + l — n) is 
large. Because a;i < 1 we can neglect the second term in 
equation (21) and {A„) cx xi". This means that the prob- 
ability of finding an A at site in the left-hand-side of the 
system is an exponentially decreasing function of the site 
index. If we write (A„) oc e~~s^ , we find A = —l/ln{xi) 
for the characteristic length of the decrease. The loga- 
rithm makes this length only a slowly varying function 
of the rate constants (see figure 2). When Wdifr becomes 
larger, a approaches 0, xi approaches 1 and A diverges. 
Note that this is a MF result. We will sec that in the 
simulations A remains finite. Also when the conversion 
is slow more A's are found away from the marginal sites. 
The second factor in the expression for a equals the re- 
ciprocal of a site being vacant. Low loading leads to a 
smaller a than high loading. Because of the vacancies 
the A's can penetrate farther into the system before be- 
ing converted. For slow conversion or fast diffusion a is 
small and A can be approximated by 



difF 



des 



The total loading with A's, Qa, is 



(32) 



(33) 
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W'ads 




W^diff 




a) 


0.2 


0.8 


0.05 


0.01 


b) 


0.2 


0.8 


0.05 


0.1 


c) 


0.2 


0.8 


2 


0.1 


d) 


0.2 


0.8 


1 


2 


e) 


0.2 


0.8 


10 


2 


f) 


0.8 


0.2 


0.05 


0.01 


g) 


0.8 


0.2 


0.05 


0.1 


h) 


0.8 


0.2 


2 


0.1 


i) 


0.8 


0.2 


1 


2 


j) 


0.8 


0.2 


10 


2 



TABLE I: The sets of parameters used for the simulations 



SO, the expression for Qa is 



n - " V^r^ n , „ s-n-«i _ 2a xijl - xi^) 



Qb = Q — Qa- 
The total production of B's is 



(35) 



Bprod = W,^QaS = 2aW, f^^^_ . (36) 



2. Simulation results 

We present now the results for different sets of param- 
eters and we compare them with MF results. Because 
we can see from equation (36) that larger pipes don't in- 
crease the productivity of the system, we consider for the 
comparisons of the results a system size S = 30. We have 
considered separately the sets of parameters in Table I. 

The sets of parameters from a) to e) are for the cases 
of low loading and from f) to j) for the high loading. The 
parameters in the table describe the following situations: 
a) and f) for very slow reaction and slow diffusion; b) and 
g) for slow reaction and slow diffusion; c) and h) for slow 
reaction and fast diffusion; d) and i) for fast reaction and 







Qa 






Q 


set 


MF 


Sim 


MF 


Sim 


Sim 


a) 


0.0330 


0.0318 


0.0099 


0.0100 


0.209 


b) 


0.0149 


0.0148 


0.0491 


0.0472 


0.198 


c) 


0.0385 


0.0342 


0.1156 


0.1024 


0.204 


d) 


0.0040 


0.0041 


0.2449 


0.2463 


0.200 


e) 


0.0046 


0.0044 


0.2767 


0.2729 


0.201 


f) 


0.0798 


0.0748 


0.0239 


0.0235 


0.795 


g) 


0.0376 


0.0373 


0.1129 


0.1157 


0.804 


h) 


0.0598 


0.0486 


0.1796 


0.1406 


0.802 


i) 


0.0048 


0.0049 


0.2931 


0.2943 


0.797 


j) 


0.0050 


0.0049 


0.3013 


0.2957 


0.801 



For the case Wads 



TABLE II: Simulation and MF results for Qa and Bprod for 
all the sets of parameters 



slow diffusion; e) and j) for fast reaction and fast diffu- 
sion. 

We can see from Table II that the simulation and MF 
results match for all the cases except the cases when we 
have low reaction rates and fast diffusion for both low 
and high loading. In these cases MF overestimates the 
amount of A's in the pipe, and consequently overesti- 
mates the B production. In figure 4 we have the site 
occupancy with A and B both from the simulations and 
MF. We again see that the MF and the simulation re- 
sults agree reasonably well, except for low reaction rates 
and fast diffusion. MF overestimates the characteristic 
length A and allows A's to penetrate farther into the 
pipe than in the simulations. The reason for this is that 
MF describes the fact that the particles cannot pass each 
other by reducing the diffusion, but this effectively does 
allow for passing. The larger A in MF means also a 
larger Qa- As a consequence the B production in MF is 
larger and, because these B's have to be able to leave the 
pipe via desorption, the probabilities (Bi) and (Bs) are 
larger in MF. The probabilities (Ai) and (A5) are there- 
fore smaller, which means that the MF curves and the 
simulation curves in figure cross each other, as can actu- 
ally be seen. The behavior of the system at high loading 
and at low loading is about the same, except that A is 
smaller at high loading. 

One might expect that the larger the number of reactive 
sites the more B's will be produced in the pipe. From the 
simulations we see that the amount of B's produced per 
unit time by all reactive sites goes to a limit value when 
the number of reactive sites is increased. In figure 5, the 
marked line represents the B production as a function 
on the length of the pipe and the dashed line the B pro- 
duction according to MF. For short pipe lengths, the B 
production from both MF and simulations increase lin- 
early with S, while for higher lengths it converges to a 
limiting value. The limiting value is higher for MF. This 
could be seen also from the Table II. According to MF 
there are more B produced in the pipe. 



B 



prod 



00 we have 

2W,xWdeB 
W„ + Wdes 



(37) 



From the simulations (see figure 3) we see that for high 
adsorption rates, Sprod converges to a point and the 
corresponding value is equal to the analytical value for 
the case adsorption is infinitely fast. The reason for 
this is that all the sites are occupied, diffusion is com- 
pletely suppressed, and only the marginal sites play a 
role. The expression above can be seen as a factor 2 
for the two marginal sites, the probability that an A at 
the marginal sites is converted to a B before it desorbs 
Wrx/VFrx + Wdes, and the rate constant for desorption 

Wdes. 

The accuracy of the simulation results for Qa and 
Bprod can be derived by looking at the total loading Q in 
Table II. For the total loading Q, the simulation results 
can be compared with the values of the exact expression 
(12). We remark that the largest deviation from the ex- 
act analytical results is 0.04, so the relative errors are 
around 0.02%. 

The differences between MF and the simulations be- 
comes especially clear in the limit Warn ^ 00. Because 
this limit makes the system homogeneous in MF we get 

^ Wads W,: 



Wads + Wdes + W, 



des 



(38) 



Wads 



Wads + Wdes Wrx + W, 



des 



(39) 



The first factor in these expressions is the probability 
that a site is occupied. The second factor indicates if the 
particle is converted to a B or not before it desorbs. The 
simulations show that the system should not be homoge- 
neous at all (see figure 6). The B production increases 
linearly with S only for the case of infinitely-fast diffu- 
sion, otherwise it converges to a limiting value. 



C. Only some of the sites reactive 

We consider now the situation that not all the sites 
are reactive, and that these reactive sites can be either 
uniformly distributed inside the pipe or distributed in 
compact blocks. We will show that the number of reac- 
tive sites doesn't change qualitatively the properties of 
the system. Qa, Qb and number of B produced for a 
variable number of reactive sites are compared with the 
previous results. 



1. Mean Field 

From the Master Equation it is easy to show that the 
total loading is again just the same as in the case when 
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FIG. 4: The site occupancy with A ((A„)) and B ((B„)) as a function on site number for cases a, b, c, d, e when S = Nreac = 30. 
The continuous line and the corresponding symmetric line represent MF results for {A„) and (B„) respectively. The dashed 
lines represent DMC results for (A„) and {B„). (A„) is decreasing towards the middle of the pipe while {B„) is increasing. 
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FIG. 5: B production as a function on the length of the pipe for Wads=0.2, WdeB=0-8, Wdis=2, and Wrx=0.1. The marked line 
represents the DMC results and the dashed line represents the MF results. 



0.9 
0.8 

0.7 




5 10 15 20 25 

n 



FIG. 6: Analytical and simulation results for site occupancy of a system when parameters are: S = Nreac = 30, Wads=0.8, 
W^des=0.2, iy(iiff=100 and Wr^=2. The continuous line and the corresponding symmetric line represent the simulation profiles 
for site occupancy with A and B particles. The bottom and the upper straight lines represent the analytical results for 
occupancy with A and with B particles respectively. 



all the sites are reactive. We introduce an extra coef- 
ficient A„ in the MF equations to the reaction term. 
An = 1 if n is a reactive site and A„ = if it is not 
a reactive site. The steady-state equations arc identical 
to equations (25), except that T^rx should be replaced 
by WrxA„. The resulting set of equations is Hnear again 
and it should be possible to solve them numerically. In 
fact only the probabilities for the marginal and reactive 
sites have to be solved numerically. For the other sites 
the probabilities can be obtained by simple Hnear inter- 
polation. That this is correct can be seen because those 
sites only have the diffusion term. We can also remove 
the probabilities for the B's, because we have from the 
model without conversion that 



(A„) + (B„) = !-(*„) 



ads 



(40) 



The resulting equations for the reactive sites have the 
same form as equation (25) for the non-marginal sites. 
We expect therefore that we get an exponential decrease 
of (A„) on the reactive sites when we move from the 
marginal sites to the center of the pipe, and a linear de- 
pendence on n between the unreactive sites. 



2. Simulation results 

The number of reactive sites is considered to vary from 
1 to 50% and the reactive sites are distributed either in 
blocks situated near the marginal sites, in the middle of 
the pipe, or homogeneously distributed in the pipe. We 
will first compare the MF results with the MC simula- 
tion results for different sets of parameters and then we 
look at the dependence of B production and total load- 
ing Qa on the number and position of reactive sites. For 
the comparison between MF and MC results we consider 
the system size S = 30 and the number of reactive sites 
Nreac = 10. The sets of parameters used for the specific 
situations to be studied are the same as the sets used in 
the case with all the sites reactive in the previous section. 

We can see from the Tables III and IV that when the 
reactive sites are homogeneously distributed or situated 
as a block in the middle of the pipe, there are significant 
differences between MF results and MC results. When 
the reactive sites form blocks near the marginal sites, the 
results are almost the same as when all sites are reactive: 
the MC and the MF results differ if we have fast diffusion 
and slow reaction. The sites in the center of the pipe are 
not relevant when the sites at the ends of the pipe are 
reactive. When the reactive sites are situated only in the 
middle of the pipe, we have deviations for all the sets of 
parameters. They are very prominent for the case when 
we have high loading, fast diffusion and fast reaction. 
MF strongly underestimates A's for all non-reactive sites, 
but we have also important deviations for high loading in 
the cases with fast diffusion-slow reaction, slow diffusion- 
fast reaction, and slow diffusion-slow reaction. This is 
happening because for high loading, the end sites will 



always be occupied by a particle A and the B's will not 
be able to get out of the pipe. In MF particles effectively 
can pass each other, so B particles are then able to get out 
of the pipe. Even for the case of low loading we still have 
deviations from MF for fast diffusion and fast reaction. 
In this case MF overestimates A's for nonreactive sites. 
For fast diffusion and slow reaction, MF underestimates 
A's for nonreactive sites and for slow diffusion and slow 
reaction MF overestimates B's for the reactive sites in 
the middle. 

Figures 7 and 8 show how the probabilities (A„) and 
(B„) vary in the pipe. The situations for reactive sites 
forming blocks at the ends of the pipe are not shown 
as they are almost the same as when all the sites are 
reactive (see figure 4). When the reactive sites are ho- 
mogeneously distributed the plots look also very similar 
to the ones with all sites reactive, except that the char- 
acteristic length A is larger. (A„) and (B„) look very 
different when the reactive sites form a block in the mid- 
dle of the pipe. The MF results show, as predicted, a 
linear behavior at the nonreactive sites. The MC results 
show, however, a nonlinear behavior in the form of S-like 
curves. At the reactive sites the behavior is similar to 
the situation with all sites reactive with the MC results 
showing a more rapid approach to the value at the mid- 
dle of the pipe than MF, i.e., smaller A. The values at 
the marginal sites can differ between MC and MF quite 
a lot. This refiects the difference in Bprod mentioned be- 
fore: a different Bprod must be accompanied by a different 
B desorption at steady-state. As we have already seen 
from the case when all the sites were reactive, Bprod very 
rapidly approaches the limiting value when the pipe is 
made longer (see figure 5). Similarly when we start with 
few reactive sites and, instead of increasing the length of 
the pipe, we increase the number of reactive sites. The 
loading Qb is already almost the same as the value with 
all sites reactive when only about 10% of all sites are re- 
active provided there are reactive sites at or very near 
the marginal sites. If the reactive sites are moved away 
from the ends of the pipe, then the loading Qb and the 
B production decreases. 



IV. SUMMARY 

We have used analytical and simulation techniques to 
study the reactivity in Single-File Systems. 

The MF results show that MF models Single-File be- 
havior by changing the diffusion rate constant, but it 
effectively does allow passing of particles. 

When all the sites are reactive, the simulation and MF 
results are very similar for all the parameters, except 
for the case when we have low reaction rates and fast 
diffusion. In these cases MF overestimates the amount of 
A's in the pipe. The amount of B produced per unit time 
by all reactive sites goes to a limit value when the number 
of reactive sites is increased. For high adsorption rates, 
Bprod converges to a point and the corresponding value 
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MF 


Sim 


MF 


Sim 


MF 


Sim 


a) 


0.0512 


0.0500 


0.0731 


0.0771 


0.0208 


0.0469 


b) 


0.0153 


0.0152 


0.0712 


0.0788 


0.0138 


0.0206 


c) 


0.0590 


0.0672 


0.0901 


0.0881 


0.0719 


0.0594 


d) 


0.0041 


0.0041 


0.0667 


0.0730 


0.0120 


0.0123 


e) 


0.0067 


0.0006 


0.0447 


0.0583 


0.0126 


0.0121 


f) 


0.0896 


0.0752 


0.3008 


0.3473 


0.1121 


0.1056 


g) 


0.0376 


0.0369 


0.2850 


0.3383 


0.0585 


0.0605 


h) 


0.0871 


0.0579 


0.2844 


0.3250 


0.1227 


0.0867 


i) 


0.0048 


0.0048 


0.2606 


0.3137 


0.0319 


0.0413 


j) 


0.0056 


0.0053 


0.1556 


0.2826 


0.0175 


0.0289 



TABLE III: Simulation and MF results for Qa for all the sets of parameters in the cases of homogeneous distribution of the 
reactive sites, blocks of reactive sites in the middle of the pipe and near the marginal sites (S = 30, Nreac = 10) 
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a) 


0.0099 


0.0099 


0.0001 


0.0000 


0.0011 


0.0045 
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0.0449 


0.0477 


0.0001 


0.0008 


0.0153 


0.0117 


c) 


0.1156 


0.1021 


0.0393 


0.0216 


0.0728 


0.0521 


d) 


0.2449 


0.2492 


0.0283 


0.0161 


0.1357 


0.0891 


e) 


0.2767 


0.2763 


0.1482 


0.0646 


0.2410 


0.1899 


f) 


0.0239 


0.0235 


0.0014 


0.0006 


0.0086 


0.0076 


g) 


0.1129 


0.1160 


0.0015 


0.0000 


0.0139 


0.1171 


h) 


0.1796 


0.1421 


0.0470 


0.0059 


0.1212 


0.0661 


i) 


0.2931 


0.2941 


0.0288 


0.0069 


0.1526 


0.0897 


j) 


0.3013 


0.2965 


0.1552 


0.0143 


0.2713 


0.1739 



TABLE IV: Simulation and MF results for Bprod for all the sets of parameters in the cases of homogeneous distribution of the 
reactive sites, blocks of reactive sites in the middle of the pipe and near the marginal sites (S = 30, Nreac = 10) 

is equal to the analytical value for the case adsorption is 
infinitely fast. The sites in the middle of the pipe have 
no effect on the B production. The differences between 
MF and the simulations becomes especially clear in the 
limit WdifF oo. 

When only some of the sites are reactive, there are sig- 
nificant differences between MF and MC results when the 
reactive sites are homogeneously distributed or situated 
as a block in the middle of the pipe. When the reactive 
sites form blocks near the marginal sites, the results are 
almost the same as when all sites are reactive: The MC 
and the MF results differ only when we have fast diffu- 
sion and slow reaction. The sites in the center of the pipe 
are not relevant when the sites at the ends of the pipe 
are reactive. When the reactive sites are situated in the 
middle of the pipe, we have deviations for all the sets of 
parameters. They are very prominent for the case when 
we have high loading, fast diffusion and fast reaction. MF 
strongly underestimates A's for all non-reactive sites, but 
we have also important deviations for high loading in the 
cases with fast diffusion-slow reaction, slow diffusion-fast 
reaction, slow diffusion-slow reaction. The MF results 
show a linear behavior at the nonreactive sites. The MC 
results show, however, a nonlinear behavior in the form 
of S-Hke curves. The loading Qb is already almost the 
same as the value with all sites reactive when only about 



10% of all sites are reactive provided there are reactive 
sites at or very near the marginal sites. If the reactive 

sites are moved away from the ends of the pipe, then the 
loading Qb and the B production decreases. 
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VI. APPENDIX 

A. Probability to find the system in a certain 
configuration. Loadings and fluctuations. 

We show the existence of a function q, depending only 
on the number of particles such that 

Pa = (41) 

is the steady-state solution of the Master Equation (2) for 
a system without conversion, where n{a) is the number 
of particles in configuration a. The second part of the 
proof consists of showing the uniqueness of the solution. 



FIG. 7: The site occupancy for the cases a, b, c, d, e -homogeneous distribution. The continuous line and the corresponding 
symmetric line represent the MF results. The other dashed lines represent the DMC results. 
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FIG. 8: The site occupancy for the cases f, g, h, i, j - middle sites reactive, Nreac = 10. The continuous line and the 
corresponding symmetric line represent the MF results. The other dashed lines represent the DMC results. 



Substitution of Pa = q{n{a)) in equation (9) shows 
that the last term in the Master Equation vanishes, be- 



(diff) 

a0 



^^Pa^^ and n{a) = n(/3). The other terms 



cause A 

can also be simplified by using how the number of parti- 
cles changes upon adsorption and desorption. 



dPg 
dt 



ads 



(ads) 



g(n(a) + l)^Ay-g(n(a))5:A(t' 

(42) 



A further simplification is possible if we realize that des- 



orption reverses the effect of an adsorption and vice versa. 
lis m 

dPa 



This means A^'';"^ = A','''''^ This leads to 

ap pa 



dt 



[q{n{a) - l)W,d. - q{n{a))Wde.] ^A 



(ads) 
a/3 



(43) 

We denote by N, the number of particles in a certain 
configuration, = n{a). We see that we get a steady- 
state solution for 



dPg 
dt 



0, 



(44) 



provided by 



qjN + l) ^ 
q{N) Wde. 



(45) 



for N = 0, 1, 2,. . . , S-1. (Note that the case A' = 5 in 
the Master Equation presents no problems, because the 
summation over j3 yields zero.) 

The second step consists of showing that this solution 
is the only one. This part for instance can be found in 
Chapter 5 of Van Kampen. [21] 



B. Derivation of function q(N) 

Expression (45) leads to 
qiN) = C 



ads 



N 



(46) 



where C is some normalization constant. We can com- 
pute it from 



a a N—O ^ ' 

S 



Wads 



Wde 



C 



q{N) 



Wdes + Wads 



Wde 



(47) 



The combinatorial factor after the third equal sign de- 
rives from the number of configurations with N particles. 
The last step uses 



The expression for q{N) now becomes 
q{N) 



Wdes 


S 


'Wads' 
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Wdes + Wads. 




[WdesJ 





(48) 



(49) 



Note that this expression does not depend on Wdiff: i.e., 

diffusion has no effect at all on steady-state properties. 

The probability p{N) that there are N particles in the 
system is given by 



p{N) 



qiN) 



W 



des 



Wdes + Wads 



W 



ads 



Wdes 



50) 

This follows from (49). With this formula we can com- 
pute all statistical properties of the number of particles. 
The average number of particles is 



E ^p(^^ 



N=0 



Wde 



Wdes + Wads 
Wads 



S S 

E 

JV=0 



N 



Wads 
Wdes 



N 



Wdes + Wads 



S. 



(51) 



The loading of the pipe, defined as the average number 
of particles per site, is 



EN^oNpjN) t^ads 



S Wads + Wde 

The average squared number of particles is 

s 



(52) 



N=0 



Wde 



s s 



E 



Wdes + Wads. 

Wads(Wdes + gWads) 
(Wdes + Wads)2 



W, 



ads 



Wdes 



JV 



s. 



(53) 



The variance, i.e., the square of the fluctuation in the 
number of particles, is then 
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(54) 



by 



C. Derivation of the one-site and two-sites 
occupancy for the model without conversion 

The probability that site n is occupied by A is given 



(A„) = ^P„A("), 

= EE^«ah, 

N aGN 
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(55) 
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where A^"^ is 1 if site n in configuration a is occupied by 
an A particle, and it is otherwise. The combinatorial 
factor denotes the number of ways the particles except 
the one at site n can be distributed over the remaining 
sites. Knowing the one-site occupancy we can derive the 
two-site occupancy 



D. Continuum limit 



(A„A„+i)=^P„AMa("+i), 

a. 

AT a&N 

= 5]g(7V)5]A(r)A(r+i), 
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(56) 



The rate equation for the A's is 

-(*„_iA„)]-H';x(A„). 



- WdifF[(A„_i*„) + (*„A„+i) - (A„*„+i) ^^^^ 



The MF approximation of this equation is 
d(A„) 



dt 



Wdiff[((A„_i) + (A„+i))(*„) 

(A„)((*„_i) + (*„+!))] 
W,x(A„). 



If we take the continuum hmit and denote by a = a{x, t), 
b = b{x,t) and v = v{x,t) the probability distribution of 
A's, B's and vacancies respectively, and if we use Taylor 
series for the diffusion term, the equation becomes 



— = Wdiff [a((n - l)d, t) + a((n + l)d, t)]v{nd, t) 



- WdifFa(nd, t)[v{{n - t) + v{{n + t)]a{nd, t) - W^^a 
« Wdiff 



dx 2 dx"^ dx 2 dx"^ 
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■ Wdiffd' 
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dx 2 dx^ dx 2 Sa:^ 



2a + d' 
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w - Wdiff a 



2w + d- 



a^a d^a d% 
aa;-^ oa;^ ax^ 



dx"^ 



a — Wrxfl 
- Wrxfl 



r 



(59) 



where d is the distance between sites. A similar relation 
can be derived for b{x, t). With D = Wdiff d-^ we can write 



J 



(da/dt\ _ n A - 6 a \ fd'^a/d. 



[db/dt ) b 1-a) [d'^b/dx'' ) + ^" 

I 



(60) 



E. MF derivation of the total loading Qa in case so, the expression for Qa is 
with conversion 

QA = ^t{A„), (61) =|E-i" + |E-i"^^-", 

= 1 E-i"' 



The total loading with A's, Qa, is written as 

s 

S 



(62) 



2a a;i(l — Xi^) 
S 1 — a;i 
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